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JULIEN GUYON 

Abstract. Given a smooth R'^-valued diffusion {Xf ,t G [0, 1]) starting at point x, we 
study how fast the Euler scheme X"'^ with time step 1/n converges in law to the random 
variable Xf. Precisely, we look for which class of test functions / the approximate 
expectation E [/ (X"'^)] converges with speed 1/n to E [/ {Xf )\. 

When / is smooth with polynomially growing derivatives or, under a uniform hypoel- 
lipticity condition for X, when / is only measurable and bounded, it is known that there 
exists a constant Cif(x) such that 

(1) E [/ (Xr-^)] - E [/ (X?)] = Cif{x)/n + O (l/n^) . 

If X is uniformly elliptic, we expand this result to the case when / is a tempered 
distribution. In such a case, E[/(Xi)] (resp. E[/(X"'^)]) has to be understood as 
{f,p{l,x,-)) (resp. {f,pn{l,x,-))) where p{t,x,-) (resp. p„(f, a;, •)) is the density of 
(resp. X"'^). In particular, ^ is valid when / is a measurable function with polynomial 
growth, a Dirac mass or any derivative of a Dirac mass. We even show that ([T]) remains 
valid when / is a measurable function with exponential growth. Actually our results are 
symmetric in the two space variables x and y of the transition density and we prove that 

dSdyPnit,x,y) - d'^d'^p{t,x,y) = d"d'^'K{t,x,y)/n + r„{t,x,y) 

for a function d^dyir and a 0(l/n^) remainder r„ which are shown to have gaussian 
tails and whose dependence on t is precised. We give applications to option pricing and 
hedging, proving numerical convergence rates for prices, deltas and gammas. 



1. Introduction and results 

Let d,r > 1 be two integers. Let {Q,J-',¥) be a probability space on which hves a 
r-dimensional Brownian motion B. We denote by J-'t = a{Bs,0 < s < t) the filtration 
generated by B. Let us give two functions 6 : R"' — > M*^ and cj : R"' — > M'^^''. We 
systematically use (column) vector and matrix notations, so that b{x) should be thought 
of as a vector of size d and cr{x) as a matrix of size dxr. We denote transposition by a star 
and define a d x d matrix- valued function by putting a = aa*. For a multiindex a E N*^, 
\a\ = ai + ■ ■ ■ + ad is its length and 5" is the differential operator • • • dx'^''-. 

Equipping with the euclidian norm ||-||, we denote by 

• ^poi(^'^) of infinitely differentiable functions / : M"^ — > M with polynomially 
growing derivatives of any order, i.e. such that for all a G N"', there exists c > 
and g € N such that for all x G M*^, 

(2) |97(x)| <c(l + ||xr), 

• C^(M'^) the set of infinitely differentiable functions / : M'^ ^ M with bounded 
derivatives of any order, i.e. such that d^f G L°°{M.'^) for all a G N'^. 
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We shall make use of the following assumptions: 

(A) For alH G {1, . . . , d} and j £ {1, . . . , r}, bi and aij belong to C^i(M'^) and have 
bounded first derivatives. 

(B) For alH G {1, . . . , d} and j G {1, . . . , r}, h and aij belong to C^iR"^). 

(C) There exists rj > such that for ah x,^ G M"', Ca{x)^ > r]\\^f. 

(C) is known as the uniform ellipticity condition. 

It is well known that, given x G M, the hypothesis (A) guarantees the existence and 
the P-almost sure uniqueness of a solution = (Af ,t > 0) of the stochastic differential 
equation (SDE) 

(3) Af = X + f h{X'f) ds + /* a{X^^) dB, . 

Jo Jo 

1.1. Motivation. Let us fix a time horizon T > 0. Without loss of generality, we can and 
do assume that T = 1. We try to estimate the law of Xf. To do so, the most natural idea 
is to approach A^ by its Euler scheme of order n > 1, say A"'^' = (X^'^ ,t > 0), defined as 
follows. We consider the regular subdivision 6„ = {0 = tg < < • • • < t^_i < = 1} of 
the interval [0, 1], i.e. = k/n, and we put Aq = x and, for all /c G {0, . . . , n — 1} and 

(4) Xr = A^f + b (aJ-) it (aJ-) {B, - B,^) . 

Then the random variable A"'^ is exactly simulatable and should be close in law of Af . 
Precisely, we measure the weak error between A"'^ and Xf by the quantities 

A?/(x)=E[/(An]-E[/(Af)] 

and we try to find the largest space of test functions / for which, for each x, there exists 
a constant Ci/(x) such that 

(5) A^/(x) = Ci/(x)/n + 0(l/n2). 

Practical interest of such an expansion has to be underlined (see, for instance, [Tj I14j). 
When dS]) holds, one can use the Euler scheme plus a Monte-Carlo method to estimate 
E [/(Af)] and then, in a time of order nA, gets an error of order l/\/iV + 1/n, where 
N stands for the number of independants copies of A"'^ generated by the Monte-Carlo 
procedure. Given a tolerance e <C 1, in order to minimize the time of calculus, one should 
then choose N = O (n^) and gets a result in a time of order . 

One can even do better using Romberg's extrapolation technique: if one runs N in- 
dependant copies (A^^"'^, A"|^^) of the couple (A^"'^, A"'^), which still requires a time of 
order riA, then computing Yld=i{'^f{-^'i^'^)~f{-^7i)) S^^^ estimate of E [/ (Af )] 
whose accuracy is of order l/^fN + l/n^, since © implies that E[2/(Af"'^) - /(Af •^)] = 
E[/(Af)] -|- 0(l/n^). Given a tolerance e <C 1, one should now choose N = O {n^) and 
gets a result in a time of order 1/e^/^. 

1.2. Previous results. Using Ito expansions, D. Talay and L. Tubaro [H] have shown 
that ^ holds when / G C^i(M'^) under condition 

(B') The biS and the cxj^'s are infinitely differentiable functions with bounded deriva- 
tives of any order > 1. 
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Hypothesis (B') is almost (B) but in (B') the functions bi and cjjj are not supposed 
bounded themselves. Using Malliavin calculus, V. Bally and D. Talay [T] have extended 
this result to the case of measurable and bounded /'s, with the extra hypothesis that X 
is uniformly hypoelliptic. If (C) holds, X^'^ and Xf have densities, say and 
p{l, X, •) respectively (in this paper, densities are always taken with respect to the Lebesgue 
measure). Then, for each pair {x,y), the authors [2] get an expansion of the error on the 
density itself of the form 

(6) p„(l,x,y) -p{l,x,y) = 7r(l, x,y)/n + 7r„(l, x, y)/n^. 

They also show that the principal error term tt and the remainder 7r„ have gaussian tails. 
Namely, they find constants ci > and C2 > such that for all n > 1 and x, y G M*^, 
|7r(l,x,y)| + |7r„(l,x,y)| < ci exp(-C2 ||x - 

Besides, V. Konakov and E. Mammen [9| have proposed an analytical approach for 
this problem based on the so-called parametrix method. If (B) and (C) hold, for each pair 
{x, y), they get an expansion of arbitrary order j of x, y) but whose terms depend on 
n: 

i-i 

(7) pn{l,x,y) -p{l,x,y) = ^ 7r„,i(l, x, ?/)/n* + O (l/n^) . 

i=l 

They also prove that the coefficients have gaussian tails, uniformly in n: for each i, they 
find constants ci > and C2 > such that for all n > 1 and x,y G M"^, |7r„^j(l, x, y)| < 
ci exp(— C2 \\x — To do so, the authors use upper bounds on the partial derivatives 

of p - which they find in [3] - and prove analogous bounds on p^s ones. 

A link with generalized Watanabe distributions on Wiener's space is exhibited in [12j . 
For the general case of Levy driven stochastic differential equations, ([5]) holds under regu- 
larity assumptions on / and integrability conditions on the Levy process, see [7l[T3]. The 
rate of convergence of the process (X"'^ — Xf ,t G [0, 1]) is given in O [6]. As for the 
simulation of densities, see for instance [8]. 

1.3. Purpose and method. Equations ([6]) and ([7|) can be seen as expansions of A"/(x) = 
E [/ (X"''^)] — E [/ (Xf )] in the special case when f = Sy, the Dirac mass at point y G M'^. 
We aim at giving a precise sense to such quantities when / is any tempered distribution, 
and at proving that expansions in powers of 1/n remain valid in this extremely general 
setting. Moreover, we will derive expansions that are valid not only for t = 1, but also 
for any time t G (0, 1], the stepsize 1/n being fixed, and we shall make explicit, in these 
expansions, the way the coefficients and the remainders depend on t, f and x. 

To get these precise results, we shall place ourselves in a strong situation. Namely, we 
will assume infinite regularity and boundedness of the coefficients of the SDE that 
is condition (B), and uniform ellipticity, that is condition (C). The reason for this is the 
following. Let us write Ptf{x) = E[/(Xf)] and = E[/(X"''')]. We first expand 

A" = — Pt as an endomorphism of C^j(M'^), in powers of 1/n. This can be done under 
nothing more than hypothesis (A) , see Theorems [9] and [10] in Section 11.81 The coefficients 
in these expansions are operators of the form Jq PgDPt-s ds or ^^n^^P^DPt-t^^, where 
-D is a differential operator. Now, under (B) and (C), both Xf and X^'^ have regular 
densities, say p{t, x, •) and Pn{t, x, •), with gaussian tails, as soon as t > 0, so that we may 
express these operators as integral operators on W^. For instance, for / G C^[(M'^), x G M"^ 
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and < s < t, 

PsDPt-sf{x) = [ pis,x,z)DPt-Jiz)dz 

= / fiy){ Pis,x,z)D{p{t - s,-,y)){z) dz] dy. 
Now the expansions read on the density itself, with coefficients of the form 
(8) vr(t,a;,y) = / / p{s,x, z)D{p{t - s, ■,y))iz) dzds, 



Jo JR'' 

(9) or 7rn{t,x,y) = V / Pn{tl,x, z)D{p{t - 4, ■,y)){z) dz. 

At this step, the key point is to prove that these coefficients, as well as any of their 
spatial derivatives, have gaussian tails (see Proposition [5]). To do so, we split the above 
time integral (resp. sum) depending on whether s (resp. t^) is small or large, and integrate 
by parts in the latter case. This is very similar to V. Bally and D. Talay's technique 
d], but they use the Malliavin calculus integration by parts formula whereas we only use 
the genuine one. Then we use upper bounds on the partial derivatives of p and as is 
done in V. Konakov and E. Mammen's work [9]. Here the uniform ellipticity hypothesis 
is crucial: it provides upper bounds that have enough quality in t to allow us to conclude. 

The same analysis, with a bit more work, can be done for the remainders. We then get 
functional expansions of the form 

i 

(10) Pn-p = T^/n + TTn/n^ Or Pn-P = '^T^n,i/n'- 

i=l 

where vr, 7r„ and the 7r„^j's and all their spatial derivatives have gaussian tails, uniformly 
in n. We then achieve to give a distributional sense to expansion ([5]) by a duality ap- 
proach: any tempered distribution can be integrated or bracketed in the variable y with 
the expansions. Theorems O [7] and [8] provide precise statements, see Section [L71 

1.4. A first series of results. Stating Theorems [6] and [8] requires a bit of preparation, 
namely defining appropriate functional spaces in which will live the coefficients vr, vr^ and 
'Kn,i in expansions (fTO|) . Before doing this, to encourage the reader, we would like to state 
a series of easy consequences of Theorem [6l including an application to financial markets. 
They will be proved in Section [1.71 The function tt which appears in them is the principal 
functional error term. It is defined by (I2ip -(|22p. Note that analogous corollaries can be 
derived from Theorem [8] as well. The first result gives the rate of convergence of the spatial 
derivatives of the density: 

Proposition 1. Under (B) and (C), for all a,/3 G N'*, there exists ci > and C2 > 
such that for all n > \, t ^ (0, 1] and x,y ^ , 

dxdyPn{t, X, y) - d^d^p{t, x,y) = ^ d°d^Tr{t, x, y) + r„(t, x, y) 

and 

|r„(t, X, y)| < cin"2t-(l"l + \P\ +'^+4)/2 f _ y||2 



The next proposition states that ([5]) is valid for measurable and polynomially growing 
/'s: 
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Proposition 2. Assume (B) and (C). Let f : -^M. be a measurable function such that 
there exists c' > and g G N such that for all x G M'*, < c'(l + Then there 

exists c > such that for all n > I, t €z (0, 1] and x G M.'^, 

(11) nfiX^)] - )] = - [ /(y)vr(t, x, y) dy + r^it, x) 
and 

\rnit,x)\ < cn^'^t^'^ (1 + ||x||^) . 

As far as extending the class of /'s for which ([5]) holds is concerned, we can even do 
better. Indeed, if for /j, £ (0, 2) we denote by £fj_ the set of all measurable functions 
/ : R'^ ^ M such that there exists ci, C2 > such that for all y G M'^, 

|/(y)| < ciexp(c2 llyD, 

we have 

Proposition 3. Under (B) and (C), for all fi G (0,2) and f G £fj,, there exists ci,C2 > 
such that for all n>\, t £ (0, 1] and x G M"^, fi^t) ^'^'^ fi-^t'^) ^''"^ integrable and 

(12) nf{Xr)] - IE[/(^f )] = - / f{y)<t, X, y) dy + r„(t, x) 
with 

l^nit, x)\ < cin~'^t~'^ exp (c2 ||a;||'') • 

In particular, ([5]) remains true under (B) and (C) when f £ £ = U^g(o,2)'?/^- More 
generally. Theorem [6] leads to 

Proposition 4. Under (B) and (C), for all a G N'^, /U G (0,2) and f G f^, there exists 
ci, C2 > such that for all n > 1, t £ (0, 1] and x G M^, 

(13) d^E[f{Xr)] - OmfiXf)] = - I f{y)dMt, X, y) dy + r^it, x) 
with 

Kit, x)\ < cin-2t-(l"l +4)/2 exp (c2 ||xf ) . 

This result can now be used in the context of financial markets. 

1.5. Application to option pricing and hedging. Let = {S^'^ , . . . , S^''^) be a 
basket of assets satisfying 

^ = ^,{S^) dt + j2 ^^AS^) dB}, 5o^'^ = > 0, 
j=i 

with (T G C^{M.'^) and a satisfying (C). Given a measurable and polynomially growing 
function (/>, we try to estimate the price Price = E[i;^(5(^)], the deltas Deltas = d'^j''E[(p{Sf)] 
and the gammas Gammajj = 9^'~''^^E[(/)(5|')] of the european option of maturity t and 
payoff (p ((ei,--- ,ed) is the canonical base of M'^). To do so, let us set x = Inv (i.e. 
x' = Inv') and Xf'' = ln(5j''*). Then X is the solution of §ii with b = fi - \\af /2 G 
C^(M'^), where ||it||^(x) = "^^j^i erf j{x). If we set exp(x) = (exp(x^), . . . , exp(x'^)) and 
f{x) = 0(exp(x)), we define a function f G £i and, since Price = E[/(Xf )], (fT2]) leads to 

Price" - Price = C7f"=X?j)/n + O (n^^t^^ ||lnv||)) , 
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where Price" stands for the approximated price E[/(X"'^)] and 

cr-Mv)= I Mu)^^^^^^^^^ du. 



>.)d Ul---Ud 



Besides, if we set Delta^ = d^^E[f{X^'^'''')] and Gamma^^- = dv''^'''E[fixf ^ 
shows that 

Delta" - Delta = C^''''^(t>{v) /n + (n'^'^/^ exp (ca ||ln u| 
Gamma" - Gamma = """""0(u)/n + O {n-'^r^ exp (ca ||ln t>| 



where 

= 1/ ^(,) ^2'^(^an.,ln.) ^^^ 
„ , ^ If , , c?n''''^^7r(t,ln't;,lnn) — l/j-,i9?7r(t, InVjlnn) 

ViVjJ(W^)d Ml •••ltd 

Eventually we have proved that applying the Euler scheme of order n to the logarithm 
of the underlying leads to approximations of the price, the deltas and the gammas which 
converge to the true price, deltas and gammas with speed 1/n, at least when the drift and 
volatility of the underlying satisfy (B) and (C), which in the context of financial markets 
seems not to be a restricting hypothesis. Note that the principal part of the error explodes 
as t tends to as t~^/'^ for the prices, for the deltas and for the gammas. 

1.6. Some functional spaces. In order to state our main results (Proposition [5] and 
Theorems [6] and E]) precisely and shortly, let us introduce some families of functional 
spaces. Functional expansions like pop will take place in such spaces. For / E Z, we first 
define QiiW'') as the set of all measurable functions vr : (0, 1] x M'^ x M*^ ^ M such that 

• for all t G (0, 1], 7r(t, •, •) is infinitely differentiable, 

• for all a,f3 G N"^, there exists two constants ci > and C2 > such that for all 
t e (0,1] and x,y £ M"^, 

(14) \d^d^TT{t,x,y)\ < cit-d^l + l^l+'^+'^/^g^pJ^.^^llx-yf /i) • 

We say that a subset B C Gi{^'^) is bounded if, in (fT^ . ci and C2 can be chosen indepen- 
dently on vr S 5. We also introduce the space ^(M'^) defined in the same way as Qi{ 
with ()14p replaced by the following two conditions: 

(15) \d^d^7r{t,x,y)\ < cit-^\''\ + \l^\+'^y^exp(^-C2\\x-yf /t"^ 

(16) (tt (t,x,x + yVi)) < cit-'^/'^exp(-C2\\y\ 



Note that we may always take the couple of constants (ci,C2) to be the same in both 
equations (fT5|) and (fT6]l . Indeed, if they hold with two couples {c[,C2) and (c",C2), they 
both hold with (01,02) if we take ci = c'^ V c" and 02 = 02 A 02- We say that a subset 
B cg{R'^)is bounded if, in ([E]) and ((I6D, ci and 02 can be chosen independently on vr S 
Note that in equation (116^ . the upper bound keeps the same quality in t, namely t""^/^, 
whatever the "number" a of times one differentiates the mapping x ^ it [t,x,x + yVt). 
This will be crucial when proving Proposition [5l 

It is convenient to extend these definitions to mappings that also depend on an inter- 
mediate time s € (0,t). To do so, let us denote by 71 the unit triangle {(s,t) G M^|0 < 
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s < t < 1} and, for I G Z, let us define Wz(M°') as tlie space of measurable functions 
p:TixW^ xW^ such that 

• for all G Ti, p{s,t, ■, •) is infinitely differentiable, 

• for all a, (3 G N'^, there exists two constants ci > and C2 > such that for all 
(s,t) e Ti and x,y E M'^, 

(17) \d^d^p{s,t,x,y)\ <cit-(l"l + l^l+'^+')/2g^pJ^_^2||x-yf a) • 

Again we say that a subset ;B C 7ii{W^) is bounded if, in (fT7|l . ci and C2 can be chosen 
independently on p G ;B. We also introduce the space 7i(R'^) which is defined in the same 
way as Hii'M.'^) with ([IT]) replaced by 

(18) \d^d^p{s,t,x,y)\ < cit-(H + l/3|+'^)/2exp(-c2||x-yf a), 



(19) \d^(^p[s,t,x,x + 



< cit '^/^exp ( -C2 \\y 



and we say that a subset B C 'H(M'*) is bounded if, in (jlSp and (jl9p . ci and C2 can be 
chosen independently on p G ;B. Again we may always choose the couple (ci, C2) to be the 
same in both equations (fT8|) and (fT9]) . Note that the upper bounds in (fTT]) . (fTSl) and (fT9|) 
are exactly the same as the ones in (fT5|) and p^ . In particular, they do not depend 
on s. 

Eventually, for 7ri,7r2 G ^(M'^), 5 G C^(M'^) and 7 G N'^, we define a function vn 7r2 
on Ti X M'^ X R'^ by putting 



(vTi 7r2) (s,t,a;,y) = / g{z)TTi{s, x, z)d] Tr2it - s, z,y) dz. 

Notation 82 means differentiation with respect to the second argument, here z. Operation 
*g^^ is a space convolution which naturally appears when developping the differential 
operator D in equations ([8]) and Q. 

1.7. Main results. We are now able to state our main results as follows. 

Proposition 5. Let Bi and B2 be two bounded subsets o/C/(IR'^), g G C^(M'^) and 'j G N'^. 
Then 

(i) {vTi *g^j 7r2|7ri G Bi,tt2 G B2} is a bounded subset of 7i\.y\(R'^) , 

(ii) {vTi o vr2|7ri G Bi,t:2 G •62} « bounded subset ofT-LiW^). 

Theorem 6. Under (B) and (C), 

(i) for all t G (0, 1] and x G M*^, Xf has a density p{t, x, •) and p G ^(M''), 

(ii) for all t G (0, 1], x G M'^ and n > 1, X"'^ /zas a density pn{t,x, •) and {pmn > 1) 
is a bounded sequence in QiW^), 

(iii) there exists it G ^i(M°') and a bounded sequence {iTn^n > 1) in ^4(1^^^) such that 
for all n > 1, 

(20) pn-p = Tr/n + TTn/n'^. 

These results are proved in Section [3.21 In Theorem [6l statement (i) is already known, 
see [4], Theorem 7, page 260, and statement (ii) has essentially been proved in |9J. As 
explained in Section [1.31 Proposition [U together with these two statements, is the key to 
derive statement (iii). 
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The function vr can be expressed in terms of p by 
(21) Ti{t,x,y) = \[ I p{s,x,z)Ll{p{t - s,-,y)){z) dzds, 



2 Jo iR<^ 

where the differential operator L2 is exphcitely given in terms of the functions a and b by 
(22) -Ll = Y, 



d. /-^ 1 \ 1 

i,j=l i,j,k=l 



ijk- 



Here, •, a^, tr, V and respectively stand for the inner product in W^, the A;-th column 
of a, the trace of a matrix, the gradient vector and the hessian matrix. In the case when 
t = I, ([2T]) agrees with V. Bally and D. Talay's expression for vr ([2], definition 2.2, 
page 100), but seems preferable because it does not involve differentiation with respect to 
t and makes explicit that the space differential operator L2 is of order less than 3, when 
V. Bally and D. Talay's operator U involves a fourth order differentiation in space. 

We shall now prove that if X is elliptic the expansion ([5]) is valid in the very general 
case when / is a tempered distribution. Let us denote by 5(M°') Schwartz's space, i.e. the 
space of infinitely differentiable functions : M"^ — > M such that x 1-^ x'^d^iplx) S L°°(M.'^) 
for all a,P S N"' (a;" stands for j;"^ •••x^'*), and let us denote by 5'(M'^) the space of 
tempered distributions. The seminorms {Nq,q £ N) are defined on S{R'^) by 



Nq{ip) = sup x^'d'^ipix) 

and the order of 5 G 5'(M'^) is the smallest integer q such that there is a c > such 
that \{S,^p)\ < cNq{ip) for ah ip £ S{R'^). Note that whenever vr G vr(t,x, •) and 

vr(i, -,!/) belong to 5(1^*^). More precisely, for B C Gi{R'^) bounded, there exists c > such 
that for ah vr E ^, t e (0, 1] and x,y £ W^, 

Nq{Tr{t, X, •)) < ct-(^+'+^)/2 ^ ||^||<?) Ng{TT{t, •, y)) < ct-(^+'+9)/2 (1 + \\yf) . 

Applying a tempered distribution 5 to ([20]) . t and x 01 t and y being fixed, we imme- 
diately deduce from Theorem [6] 

Theorem 7. Under (B) and (C), for all S G 5'(M'^), there exists c > such that for all 
n>l,t£ (0,1] and x,y £ W^, 

{S,pn{t,x,-)) - {S,p{t,x,-)) = -{S,TT{t,x,-)) +r'^{t,x), 

n 

{S,Pn{t,-,y)) - {S,p{t,-,y)) = -{S,TT{t,-,y)) +r'^{t,y), 

n 

and 

|r;(t,x)| + |<(t,x)| < cn-V('^+4+#^)/2 ^ ^ ||^||#5^ _ 

Let us define E [^(y)] by {S,py) when S £ 5'(M'^) and y is a random variable with den- 
sity pY £S{R'^). Note that, when S is a measurable and polynomially growing function, 
this definition coincides with the usual expectation. We then have proved that, under (B) 
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and (C), ([5]) is valid for /'s being only tempered distributions, and not only for t = 1, but 
also for any time t G (0,1], and we have even precised the way the 0(l/n^) remainder 
depends on t, f and x. Precisely, this remainder grows slower than as x tends to 

infinity, and explodes slower than ag t tends to 0. 

We can now prove the propositions stated in Section II. 4i Proposition [1] is immediate 
from Theorem |6l In the special case when S" is a measurable and polynomially growing 
function, we get Proposition [2) 

Proof of Proposition Multiplying (j20p by f{y) and integrating in y leads to (jlip with the 
remainder r„(t,x) = n~'^ f^a f{y)iTn{t,x,y) dy. Since \f{y)\ < c'(l + ||y||'^) and (7r„,n > 1) 
is bounded in ^4(R'^), we can find ci > and C2 > such that for all n > 1, t G (0, 1] 
and X G R"^, |r„(t,x)| < cin-H-'^'^+^'>/^ J^a{l + WvW^) exp{-C2\\x - yf /t) dy. Setting 

C = {y-x)/Vi leads to |r„(t,x)| < an-H-^ f^,{l + \\x + (ViW^) exp{-C2\\Cf) dC- To 
complete the proof, it remains to observe that there exists c > such that for all t £ (0, 1] 
andx,CGM'^, + C\/i|r < c(||a;f + ||Cf ). □ 

It is easy to adapt the preceding proof to get Proposition [3l In the same way, differ- 
entiating dS]) a times in x, multiplying by f{y) and integrating in y leads to Proposition 

H 

Expansion (j2U|) should be seen as an improvement of ([U]): it allows for infinite differen- 
tiation in X and y and also precises the way the coefficients explode when t tends to 0. 
We have an analogous improvement for expansion ([7]): 

Theorem 8. Under (B) and (C), for each « > 1, there exists a hounded family (vr^^j, n > 1) 
in ^2i-2(R'^) and two hounded families (vr^ j,n > 1) and (vr^ j,n > 1) in ^2i(R'^) such that 
for all j,n > 1, 

(23) Pn-P=T.^-^ + t(^ 

1=1 1=2 

Here and in all the sequel we use the convention that a sum over an empty set is zero, 
and [nt\ denotes the greatest integer less than or equal to nt. Expressions involving 
[nt\ do not appear in (j20p since they are hidden in the remainder. When t = 1 and no 
differentiation is applied neither in x nor in y, (j23p boils down to the result of V. KONAKOV 
and E. Mammen [9]. Again note that (j23p is much richer in the sense that it allows for 
infinite differentiation in space and also precises the dependence on t. Theorem [8] will also 
be proved in Section [3.21 

1.8. A preliminary result. As explained in section [L3l in order to prove point (iii) in 
Theorem [6l we first seek an expansion for the error operator 

Ar = pr - Pt 

where, for / G C;;^i(M'^) and x G R'^, we have set Ptf{x) = E[f{Xf)] and PPfix) = 
E[/(X"'^)]. Precisely, we look for operators Ct and such that i?" = 0(l/n^) and 
A" = Ct/n + R^. The following theorem, interesting in itself, is proved in Section [2j It 
can be seen as an improvement of |14) . It not only gives explicit formulas for Ctf{x) and 
R'^f{x) but also provides useful information about their dependencies on n,t,f and x. 
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Note that it does not require neither (B) nor (B') nor (C). In order to state it shortly, let 
us 

• denote by £ (C~i(M'^)) the space of endomorphisms of C^j(M ), 

• say that a subset B C C^j(M'^) is bounded if, in (l2|), c and q can be chosen 
independently on f G B, 

• say that T G £ {C^^(R'^)) is bounded if for all bounded B C Cl^ilM'^), {Tf\f £ B} 
is a bounded subset of C^^^ilM'^), 

• denote by £b {C^ii^'^)) the space of bounded endomorphisms of C^j(M'^), 

• say that a Cb (C^i(M'^)) -valued family (Ti,i G I) is bounded if for all bounded 
B C C;;^i(M'='), {Tiflf GB,ieI} is a bounded subset of C;;^i(M'^), 

• say that {Ti,i G /) is a 0(/i(i)) family in Cb {C^ii^'^)) if the family {hii)-^Ti, i G /) 
is bounded. 

It is already known that, under (A), {Pt,t G [0, 1]) is a bounded family in Cb (C~i(M'^)) . A 
proof can de found in [llj . Lemma 3.9, page 15. Using Lemma \25\ this proof straightfor- 
wardly adapts uniformly in n so that (Pj", t G [0, 1], n > 1) is also bounded in Cb {C^ii^'^)) ■ 
We are now in the position to state the main result of the first step: 

Theorem 9. Under (A), (Af,t G [0,1], n >l) is a 0{t/n) family in Cb {C^i{R'^)) , and 
there exists a 0{t) process {Ct,t G [0,1]) and a 0(l/n^) family {Rt ,t G [0,1], n > 1) in 
Cb (C;^i(M'^)) such that 

Moreover, Ct is explicitely given in terms of {Pt,t G [0, 1]) and of L*^ (see [2^l ) by 

(24) Ct = l f PsLlPt^.ds. 

^ Jo 

Note that this theorem covers the result of D. Talay and L. Tubaro [14j since it 
imphes that for any / G C|^i(M°') we can find a g G N such that for ah x G M'^, t G [0, 1] 
and n > 1, 

Ar/(x) = a/(x)/n + ofi^'' 

It even improves it a bit since we see that this holds under nothing more than condition 
(A), whereas D. Talay and L. Tubaro state their result under the stronger condition 
(B'). Note also that if we restrict ourselves to times t belonging to the discretization grid 
6„, we get a better control, of order 0{t/n'^), of the remainder, see Remark 1141 
Instead of Theorem [9l in order to derive Theorem [8l we shall need 

Theorem 10. Under (A), there exists a sequence of differential operators (L'jjj > 2), 
recursively defined by [27\)- f2^) . and for each i > 1 a 0{t/n^) family i?* = G 
[0, 1], n > 1) in Cb {C^iO^'^)) such that for all j > I, 

1=2 ■ fe=0 1=2 

Observe that the main term in (I25D is 



n \ In ^-^ ^ I n 

k=o 
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and the remainder is of order 1/v? . Note also that if we restrict ourselves to times belonging 
to the discretization grid S„, we get the following expansion in Ch (C^i(M'^)): 

i ^ L«*J-i 

^fntj/n = ^ PtlL*P[nt\/n~tl + 

i=2 ' k=0 

Theorem [TOl is also proved in Section [2j 

1.9. Organization of the paper. Section [2] deals with the expansion for the expectation: 
it is dedicated to the proofs of Theorems and [TDJ 

Section [3] is our second and final step. It is devoted to the proofs of Theorems [6] and El 
It begins with the proof of Proposition concerning the space convolution *g^^ in Q{W^). 

Eventually, Section [His an appendix where we have gathered useful results on the Euler 
scheme and technical lemmas that are used in Sections [2] and El 

2. First step: expansion for E[f{X^'^)] 

In this section we seek to expand A'^f{x) = E [/ (X"'^)] — E [/ (Xf )] in powers of the 
time step 1/n when / is a regular function, say / € C^i(M'^). The idea is the following. 
Recall the discussion preceding Theorem [9l under (A), both Pt and are endomorphisms 
of C;;^i(M'^). In C (C;^i(M'^)) we then write 

[nt\-l 

(26) Ar = Pr-Pt= Ptl^l/nPt-tl,. + Pvnt\/n^tvnn/n- 

k=0 

There is a subtle point here: {X^'^,t £ [0, 1]) is not a Markov process, since the future 
of X"'^ depends on the past value X'^^^^^, see (jl]). Nevertheless, it is easy to check by 
conditioning on J^t" that we have PJliP^ = PJ^+s all s > - but beware: this is different 
from PgPPi as soon as ns is not an integer. 

Equation (j26p leads us to expand A" for small t, namely for t < 1/n. This naturally 
involves a series of differential operators as we shall now see. 

2.1. Operators associated with the Euler scheme. Let us denote by L the infinites- 
imal generator of the diffusion X and by {L^,x G M'') its tangent infinitesimal generator, 
i.e. 

d d d 1 

L = Y, b^^'^ + 2 E ^i^J^^'^"' and = ^ bi{x)d'' + - ^ a,,, (x)r»+^^ . 

i=l *ii=l *=1 id=^ 

We use the convention that L and act on the variable y, so that, for instance, Lij){t, x, y) 
and L^ip(t,x,y) respectively stand for L {tp(t, x, ■)) (y) and {il){t,x.,-)) {y). is the 
infinitesimal generator of the Euler scheme {X^'^ ,t £ [0, 1/n]) starting from x, over the 
first discretization time interval: is built from L in the same way as X" is built from X, 
by freezing the drift b and the volatility a to their initial value on discretization intervals. 
Besides, for each x G M'^ we define a sequence of differential operators G N) by 

putting Lq = I (the identity operator) and 

(27) LJ+i = L^LJ - L^L, 
and we set 




(28) 



L*f{x) = Ly{x). 



12 



JULIEN GUYON 



Observe that L\ = 0. Besides, is given by ([22]) so that, under (A), G (C;^i(M'^)) 
and there exists a family 1 < |a| < 3) in C|^j(M'^) such that 

(29) ^2= E 52,a9". 

l<|a|<3 

gives the exact principal error term in the expansion of Af , see (|24p and (j2ip . L* is the 
differential operator appearing in ([25]) . It does not give the exact expansion in powers of 
1/n but an approximated version, in the spirit of [9], since in (|25l) the coefficients depend 
on n - but should themselves be expanded in powers of 1/n. See [7J, equations (6.35) and 
(6.36), for an expression of the operators involved in the exact expansion. 

Under (A), L and belong to £{, (C'^i(M'^)) for each x G W^, and, by induction, so 
does LJ. We can describe L^' more precisely. Indeed, defining the powers of an operator 

^ by ^0 = I and Ai+^ = AA^ , inductions on j lead to = Ei=o(-l)*(i) {L'^Y'' and 
to the existence of a family {gj^a, hj^a,j G N*, 1 < |a| < 2j) in C^i(M'^) such that 

Vx G {ly = 9jA^)d'^ and U = ^ /ij>5". 

l<|a|<2j l<|a|<2i 

Hence, for each j € N* one can find a family {mj^a, 1 < |o| < 2j) of integers and a family 
i9j,a,l, hj,a,l, 1 < |a| < 2j, 1 < / < rrijA in C^iiW^) such that for all x G W^, 

(30) q= Yl E5..",K^)/ijWp"- 

l<|a|<2j V 1=1 J 

Remark 11. Note that when (B) holds, the functions gj,a,h hj^a,i and all belong to 
C^(M°') (in fact they are polynomial in b, a and their derivatives). 

We are now in the position to define a family of operators = {^^'/,n > 1,0 < s < 
t < 1/n) as follows: 

(31) yfGC^^M''), <^:;//(x) = E[LjPi_,/(Xr)]. 
Observe that ^g'/ = L*Pt and that, from (f30]) . 

'^j,a. 

(32) Y Y9J,c.,iP^ihc.,id''Pt^s). 

l<|a|<2j 1=1 

Boundedness is a key property of this family: 

Proposition 12. Under (A), is a bounded family in Cb {C^i{R'^)) . 

Proof. {Pt,t e [0, 1]) and {Pp,t £ [0, l],n > 1) are bounded families in Cb (C~i(M'^)), see 
the discussion preceding Theorem[9l Besides, multiplication by a function in C^i(M'^) and 
differentiation are bounded operators on C|^[(M'^). As a sum of compositions of bounded 
famihes in Cb {C^i{R'^)) , is a bounded family in Cb {C^i{R'^)) . □ 

The family naturally appears when we recusively use Ito's formula to expand A" for 
small t, as we now explain. 
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2.2. Ito expansions. We recall (see [H], theorem 3.11, page 16) that for / G C^i(M'^), 
{s,y) 1-^ Pt-sfiy) is infinitely difFerentiable on [0,t] x M.'^ and 

(33) V(s, y) G [0, t] x (9, + L)Pt^sf{y) = 0. 
Since 0^ and commute, and the definition of L| imply 

(34) (5, + L-)L-jPt^s = {L-L-j - L]L)Pt^, = L]^,Pt^s- 

For a measurable family (As) in (C^j(M'^)), we denote by f^^ Ag ds the element of 
C {C^^iiR"^)) which maps / to X i — > J"^ ^ ^s/(x) ds. The following lemma states that '^. '/^ 
is the derivative of on the interval [0,t]. 

Lemma 13. Under (A), for all j G N, n > 1 and < s < t < l/n, 



(35) 



Jo 



Proof. For / G C;^i(M'^), {s,y) ^ L^Ptsfiy) is infinitely difFerentiable on [0,t] x so 
that we can apply Ito's formula to it and to the semimartingale X"'^ between and s. 
Using for the second equality, we get 

L^Pt-J (Xf '-) - L^Ptf (x) - Ms 

= f + L^) ^Pt-s'f {Xy) ds' = £ ^^,P^-,f [Xy) ds' 

where M, = ^ti E,'=i ^iA^) Jo (^^jPt-s'f dB^,. Since {LJPt_,,/|s' G [0, t]} 



is bounded in C^i(M'^), (I59p implies that (Ms,s G [0,t]) is a square-integrable martingale 
and thus has zero mean. Hence, taking expectations and using (13ip and Fubini's theorem, 
we have 

C'//(x) - L*PJ(x) = rE[L^^,P,_,>f{xy)] ds'= [\:,f'f{x)ds', 

Jo Jo 

which concludes the proof. □ 
For t G [0, l/n], since = ^>"f - by iterating ([35]) we get 

(36) Ar = ^^L*P, + /r'^+\ 



i=2 



where 



(37) ly^' = f r... r $^.^+1 ds,^,...ds2ds,. 

Jo Jo Jo 

The crucial point here is that, by construction, = so that the sum in ([36l) begins with 
i = 2. 

Injecting this in (j26p . we eventually get for all t G [0, 1] and n > 1 

j . [nti-l j . I , I / y 

(38) ^t=E^T. Pk^P^-^ + + E ^ ^n^J/n^r^.-M/n. 

i=2 ' k=0 i=2 
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where 

[nt\ - 1 

(39) i^r^' = E P^yi^'P^-^".,. + Plntl/ni:tl\/n- 

k=0 

Prom Proposition [El > l,t S [0,1/n]) is a 0(t^+i) family in Cb {C^^iR'^)). 

Recalling the boundedness of {Pt,t G [0,1]) and G [0,1], n > 1), we get that the 

family R' = {Rt'\t e [0, l],n > 1) is 0{t/n') in Cb {C^ii^'^)). Theorem[TO]is thus proved. 
We are now in good position to prove Theorem [9l 

2.3. Proof of Theorem [9l In the particular case when j = 1, ()38p reads A" = i?"'^ so 
that we have proved that {Af,te [0, l],n > 1) is 0{t/n) in Cb (C~ (M'^)), which was the 
first statement of Theorem [9l 

In the particular case when j = 2, if we set 

1 /■* 

(40) Ct = - PsLlPt^s ds, 



2 



[ntj-l 



(41) Ai, = ^ i E p^i^'^Pt-ti - 1 Psm-sds 

\ k=0 

^ M-1 

(42) Alt = ^ Yl {Pq^lPt-tl - Ptl^lPt-tl 

k=0 

f,o\ An _ jyn,2 , jt - [ntj/n)^ 

(43) ^3,t - H ^ Plnt\/n^2Pt-lnt\/n 

(44) = Al, + Al, + Al„ 
equation (138]) reads 

(45) A^ = Ct/n + R^. 

As a composition of bounded families, {PsL2Pt~s,0 < s < t < 1) is a bounded family in 
Cb {C^^{R'^)), so that {Ct,t E [0, 1]) is 0{t) in Cb {C^^{R'^)). It remains to prove that 
{R]^,t G [0, l],n > 1) is 0(l/n2) in A (C;^i(M'^)). We have already proved that it is true 
of t G [0, 1], n > 1). It is obviously also true of {{t - [nt\ /nf Pp^^j/„L^Pf_L„<j/„, t G 

[0,1], n > 1), so that {Alt,t€ [0,1], n > 1) is ©(l/n^) in Cb {C^i{R'^)). 

For {Alt,te [0, l],n > 1), observe that, if we set L* = LL^ - L^L G Cb {C'^^^{R'^)) , as 

dsPs = LPs = PsL, we have dsPsLlPts = PsL^Pts - PsL%LPt-s = PsL*Pt-s- Hence 
the family {Pt^^LlPt-ti - PsLlPt-s,t G [0, 1], n > 1, A; G {0, . . . , [nt\ - l},s G [t^,t^+i]) 
satisfies 

(46) PtlLlPt-ti - PsLlPt-s = - r PuL*Pt-u du 

and thus is 0(l/n) in Cb {C'^^^{R'^)) . As a consequence, 

1 L«*J^i ftl^, ^ , 1 /■* 

(47) Alt = ^ {P^L^Pt-ti - PsLlPt-s) ds-— PsLlPt-s ds 

fc=Q Jtl J\nt\/n 

isO(l/n2) in Cb (Cp-(M'^)). 
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As for {Alt,t e [0,1], n > 1), note that Plh^Pt^q - Pt^L^Pt-q = A^^nLlPt-q. 
Since (AJ^,t G [0,1], n > 1) is 0(l/n) in Cb {C^iiM"^)) , so is the family [P^nLlPt-ti - 
PfnL2Pt-t^,t £ [0,1], n > l,k £ {0,...,[nt\ — 1}), as the composition of a bounded 
family by a 0{l/n) family in £f, (C^i(M'^)). This completes the proof of Theorem [9l 

Remark 14. It is noteworthy that the family {R't^,t G [0, l],n > 1) defined by 



J L^^J / ^ 

is 0(t/n2) in Cb {C^^y{W^)) . In particular, G [0,1], n > 1) is 0{t/n'^) in 

Cb (Cp-(M'^)). 

3. Second step: expansion for the density of X^'^ 

This section is devoted to the proofs of Theorems [6] and [HI 

3.1. Space convolutions. We begin by proving Proposition [5] which is the key argument. 
Recall the definitions of Section [L6l Let Bi and B2 be two bounded subsets of Q{W^)^ 
g £ C^{R'^) and 7 E N°'. We want to prove that 

(i) {vri *g,^TT2\TTl £ Bl,TT2 £ B2} is a bounded subset of 7Y|^|(M ), 
(ii) {vTi *g '^2^1 G ;Bi, G B2} is a bounded subset of Tl{M.'^). 
The functions vri *g^^ tt2 depend on {s,t,x,y). We shall proceed differently depending 
on s is small or large with respect to t. The main trick is to integrate by parts in the latter 
case, so that the derivatives should always rest on the regularizing part of the integral. 
This is analogous to V. Bally and D. Talay's use of Malliavin calculus integration by 
parts formula [Ij. This is the reason why we partition the unit triangle Ti into = 
{{s,t) £ Ti|0 < s < t/2} and T+ = {(s,t) £ Ti\t/2 < s < t}, and, for e = ±, we 
define (tti *g^-y 712)^ (s, t, x, y) = lT{-is, t) (tti *g^'y 1^2) (s, t, x, y). We then have vri *g^-y tt2 = 

(tTI *g^-y 7r2)_ + (VTl *g^^ 1^2) + - 

Before proving Proposition [5] and for the sake of clarity, let us state apart the following 
technical lemma, whose proof is a straightforward application of Lebesgue's dominated 
convergence theorem: 

Lemma 15. Let I £ Z, (xi, ^ £ I) be a family of measurable functions mapping 7i x x 
M"' X M"' into M such that 

• for all i £ I, {s, t) £ Ti and C £ W^, Xii^i "i "> C) infinitely differentiable, 

• for all a, /3 £ N'^, there exists two constants ci > and C2 > such that for all 
i £ I, (s, t) £ Ti and x,y,C £ R'^, 

(48) \d^d^Xi{s, t, X, y,C)\< + l^'l exp (-C2 \\x - yf /t - C2 ||C| 

and let us define I {xi)is,t, x,y) = f^dXi{s,t,x,y,C) dC. Then {Z{xi)\i £ 1} is a bounded 
subset ofHiiW^). 

Proof of Proposition\^(i). It is enough to show that both B^ = {(tti vr2)e|7ri £ Bi,tt2 £ 
B2} are bounded. 

Step 1. Let us first treat B-^, i.e. the case when s is small. After the change of variables 

z = x + C^/s, we get (vri *g,j t^2)- = ^^(x^^i.Tra) with 

X,^i,7r2('5>i>a;>y)C) = Irf (•5>i)'5'^^^5'(a:^ + CVs)Tri{s,x,x + C\fs)dl'K2{t - s,x + CVs,y). 
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It is enough to check that the family (^Xni 712' (^11^2) & Bi x B2) satisfies the assumptions 
of Lemma[T5]with / = |7|. The first point is obvious. In order to check the second one, let 
us fix a, /9 G N'^. According to Leibniz's formula, d^dyXm, 772(8, t, x, y, C) can be written as 
a weighted sum of terms of the form 

9^ (vri(s, x,x + CV^)) d]+'''d^7T2{t -s,x + CV^, y), 



with |q!i| + \a2\ + lets] = |a|, so that in order to check (j48p it is enough to show that for 
each such (01,02,0^) one can find ci > and C2 > such that for all (vri,7r2) £81x82, 
{s,t) £ Ti and x,y,C £ M'^, 2'"^'"^(s, t, x, y, C)| is less than the r.h.s. of (gS]), with 

/ = I7I. Now, 81 and 82 are bounded subsets of Q{W^) so that from (fT5]) - (fT6]) one can find 
C3, C5 > and C4 > such that for all (vri, 7r2) £ 81 x 82, {s, t) £ Ti and x,yX £ 



\d^^ {tTi{s,X,X + CV^)) \ <C3S-^/2 

exp ( — C4 ll'^ll 



and 

V- {s, t) l^r^^Sf vr2(t -s,x + CV~s, y) 

< 1^- (s, t)c3{t - s)-(l°3l + 1^1 + 1^1 exp (-C4 \\x-y + C^||^ /{t - s)) 

< V- (s, t)c5r + 1^1 + 1^1 ||x - 2/ + C^sf /t) 

where, for the last inequality, we have used the fact that when (s,t) £ , t/2 < t — s < 
t < 1. Now, using the fact that ||x — z\\'^ > ||rE||^/2 — ||z||^ for all x,z £ M*^, we see that 
for all (s, t) £ Tf , ||Cf + \\x-y + CV^f /t > (||x - yf /t + ||Cf )/2. Since g £ C^{R''), 
we can eventually find ci > and C2 > such that for all (vri,7r2) £81 x 82, is,t) £ 71 
and x,y,C ^ 

\Xni%T'''"' ^' ^' 2/' C) I < + I'^l +'^+1^1)/^ exp (-C2 ||x - yf /t - C2 ||Cf ) , 

which completes Step 1. 

Step 2. Let us now treat 8-^, i.e. the case when s is large. After I7I integrations by 
parts, we have 

{Tri*g^^Tr2)+{s,t,x,y) = l^+{s,t) / d2{g{z)TTi{s,x, z))7r2{t - s, z,y) dz. 

Using Leibniz's formula and making the change of variables z = y — Q^Jt — s, we get that 
(tti *g^^/ 7^2)+ is a weighted sum of terms of the form 1{xtun2'^^) with 

X+'7^^f^(s,t,x,y,C) = lr+is,t){t - sf/^d^^g{y - CVt^s) 

dfT^i{s,x,y- CVt - s)iT2{t - s,y - (Vt - s,y) 

and I71I + I72I = I7I, so that we are now in the position to apply the same arguments as 
in Step 1 and get that the family {xt^nC^ , ('?ri,7r2) £ 8iX 82) satisfies the assumptions of 
Lemma [T5] with / = I7I, which completes the proof. □ 

Proof of Proposition\^(ii). From (i), we know that {vri *g,o 7r2Ki £ 81,112 £ 82} is a 
bounded subset of TioiW^). It remains to prove that (jl9p holds for p = tti *gfi 7^2 with 
constants ci and C2 which do not depend on (7ri,7r2) £ 81 x 82- As in the proof of 
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Proposition [5]- (i), we treat (vri *gfi7r2)^ and (vri *3,o^2)+ separately but analogously. That 
is, after integrating by parts, the term (vri *g,o 7^2)+ can be treated in the same way as 
{Tri*gfi'7T2)~- Thus we shall only deal with the latter term. We have {Tri*gflTr2)- = I{Xni -^2) 
with 

Then we write 5° (x^^ ,^2 t,x,x + yVi, C)) as a weighted sum of terms of the form 

x;;",™(s,t,x,y,C) = lr-{s,t)s''/^d'^'g{x + CV^) 

(vri(s,x,x + C\/s)) c?"^ {^2{t - s,x + C\/s, 2; + yVt)^ , 

with + \a2\ + jasl = Then we use (jl6p twice and the same arguments as in the 
preceding proof to get ci > and C2 > such that for all (vri,7r2) £ Bi x B2, {s,t) £ Ti 
and x,y,C & ^'^1 IXtti^'^^'^^I-s, t, x, y,C)\ ^ cit"'^/^ exp(— C2 ||y||^ — C2 ||ClP)) and an obvious 
adaptation of Lemma [15] completes the proof. □ 

3.2. Proof of Theorems [6] and [8l In this section, we assume (B) and (C). We first 
want to prove Theorem [6l We recall that statement (i) is already known, see [3], theorem 
7, page 260. The next lemma is statement (ii). 

Lemma 16. Under (B) and (C), for all t G (0, 1], n > 1 and x E , X^''^ has a density 
Pn{t,x, •) and {pn,n > I) is a bounded sequence in Q{W^). 

Proof. It is known that for all n > 1, /c G {1, . . . ,n} and x G R"', X^"'^ has a density 

k 

Pn,k{x, •) such that p^^k is infinitely differentiable and satisfies (fT5|) -(fT6 |) with t = and 
two constants ci and C2 which do not depend on n and k (see the proof of theorem 1.1, page 
278, in [9]). Since [nt\/n > t/2 for all t > 1/n, this shows that the sequence {pn,n > 1) 
defined by Pn{t, x, y) = '^{nt>i}Pn,\nt\ {x, y) is bounded in Q{W^). If we denote by T{t, x, ■) 
the density of x + b{x)t + a{x)Bt {t G (0, 1]), we observe that when k G {1, . . . ,n — 1} 
and i G (tfc,tfc+i), X"'"" has the density p„(t, x, •) = j^aPn,k{x,z)T{t-fl,z,-)dz = (p„*i,o 
r)(t^,i,x, •)• Hence, for ah t G (0, 1], n > 1 and x G W^, X"'^ has the density 

{Pn,k{x,-) if t = t^, A; G {1, . . . ,n}, 

T{t,x,-) if tG(0,ty), 

(Pn*i,or)(t^,t,x,-) if t G (t^,t^+i), /cG {l,...,n-l}. 

Observing that T G ^(M'^) and applying Proposition [5]-(ii), we get that {pmn > 1) is a 
bounded sequence in ^(M*^). □ 

We shall now prove statement (iii) of Theorem [6l Recall (j45p . We want to make explicit 
Ct and i?" as integral operators on M°'. To this end, note that, applying recursively 
Lebesgue's dominated convergence theorem, we have that for all t G (0, 1], / G C;^i(M'^), 
X G M'^ and a G N'^, 

(49) d^Ptf[x)=f fiy)dMt,x,y)dy. 

The next lemma explicits Ct as an integral operator. The function vr which appears 
there should be thought of as the kernel of C. 
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Lemma 17. Under (B) and (C), there exists vr G Qi{W^), given by [21\) . such that for all 
t£{0,l], f £ C^.iiR'^) and x G M^, 

Ctf{x)= f{y)TT{t,x,y)dy. 



Proof. Using (j40p for the first equality, (j29p for the third one and (|49p for the fourth one, 
we have 

2Ctf{x) = f PsL*2Pt-sf{x)ds 
Jo 

/ p{s, X, z)L*2Pt^J{z) dzds 
rt 



y2 f f 9l,a{z)p{s,x,z)d'^Pt-sf{z)dzds 
y2 f f f fiy)92,a(.z)p(.s,x,z)d2p{t-s,z,y)dydzds. 



l<|a|<3 



Using Fubini's theorem, we see that to complete the proof it is enough to show that the 
function vr defined by 

(50) 7r(t,x,y) = - ^ / {p *g*^^^a p){s,t, x,y) ds 

l<|a|<3 ° 



belongs to Now, p G g{R'^) and, from Remark [IH g^^^ G C^°°(M'^) so that we can 

apply Proposition [5]-(i): p*g*^,aP G T~(-\a\(R'^) ■ In particular, Jq{p *g* p){s, ■, ') ds G 
^|o|-2(^°')- Since \a\ < 3 and by monotonicity of {Qi{W^),l G Z), we finally get that 
vr G t/i(M'^). To complete the proof, note that (jSOp can be rewritten as (j2ip . □ 



We have a similar representation for A^^, recall (j4ip . We say that a sequence (vr"", n > 1) 
is 0(l/nJ) in Gii'R'^) if (n%",n > 1) is bounded in ^^I^'^)- 

Lemma 18. Under (B) and (C), there exists a 0{\/n^) sequence (vr",n > 1) in ^3(M'^) 
such that for all t G (0, 1], / G C;;^i(M'=') and x G , 

Altf{x)= [ f{y)n^{t,x,y)dy. 

Proof. Recall ()46p . From Remark II H there is a family {gf^, 1 < |a| < 4) in C^{R'^) such 
that uf = ^i<|a|<4 ^^qC^", so that, using ([19]) . we have 



{PtnLlPt-ti-PsLlPt^s)f{^ 



l<|a|<4 



fiy)93ai^)pi^^ 2;, z)d2p{t - u, z, y) dydzdu. 
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Using (HZ]), we get AlJ{x) = f^^ fiv^iit^ ^> v) dy with tt^^ = 7rf_i + t^^^ and 
Now Proposition [5]-(i) states that p ^ p and p *g*^,a P belong to Tl\a\{^'^)- Hence 

in+' ffnip** 
-'^k 93,0 



ilt^^^ It^iP* * •,-,•) duds,n > 1,/c e {0,... ,n - 1}) is 0(1^) in a|„|(M'') and 



(Iln-\/niP*92,c'°'P^^^''''''^ds,n > 1) isO(l/n) in ^|q,|(M''). As a consequence, {iT'^^^,n> 1) 
is 0(l/n2) in ^2(1^'^) and «_2>™ > 1) is 0(1/^^) in ^3(M^). Eventually, (7r^,n > 1) is 
0(l/n2) in g3(M'*). ' □ 

We shall now prove an analogous lemma for ^2 t- 



Lemma 19. Under (B) and (C), there exists a 0(l/n^) sequence {1^2, n > 1) in ^3(M° 

~fOO I 

-'pon 



such that for all te{0,l], / G C^oii^"^) and x E M^, 



Altfi^)= / f{y)7r^{t,x,y)dy. 
Proof. Since P^^LlPt^t^^ = PqL*2Pt^q when /c = 0, (gl]) reads 

[nt\ - 1 
A:=l 

[Jlt J - 1 



E / {Pn-pm,X,z)LlPt^t-f{z)dz 

E E / iPn-pM,x,z)glJz)d''Pt-t^J{z)dz 

:\a\<3 k=l -^^^ 
[nt\~l 

EE// {Pn-p){tl,x,z)gl^{z)f{y)d^p{t-tl,z,y)dydz 



l<|a|<3 k=l 



where we have used (|29p for the third equality and ()49p for the fourth one. From Remark 
dH G Cft°°(M'^) so that to complete the proof it is enough to show that whenever 
g G C^{R'^) and a G N"', the sequence (7r",n > 1) defined by 

7r"(f,a;,y)= V / iPn-p){tl,x, z)g{z)d^p{t-tl, z,y)dz = V {{pn-p)*g,aP){tl,t,x,y) 

is bounded in ^|q,|(M^). And to do so, it is enough to show that the sequence (/o"n,n > 
2, /c G {1, . . . , n - 1}) defined by 

p%{t,x,y) = lri(ifc,i) {{Pn-p) *g,aP) itk,t,x,y) 
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isO(l/n) in _,_2(I^'^)- Let us write /3"n~(t, x, y) = l^~{t'^,t)pfn{t,x,y) and p^,l^ {t, x,y) = 

' ' k Ik k 

\+itk^t)pfn{t,x,y) so that = p^f + Ptj^ ■ 

Let us first prove that {p?A~,n > 2,k e — 1}) is 0(l/n) in G\a\+2(^'^)- The 

k 

sequence (vrtn,n > 2,k G {l,...,n — 1}) defined by Trt^{t,x,y) = 1^- {t^,t)g{x)d2p{t — 
t^,x,y) is bounded in ^|q,|(M'^), since t — > t/2 when G T^^. Now note that 

p^,i~ = P^T^ti — Pt^T^ti = '^i^'^ifc (s^^ ([^T|) in the appendix for the definition of P"vr, P^tt 
and A> when vr G ^Kl^'^))- Thus, from ([MD-dSSD and ([37]) apphed with j = 1, 

Proposition [23] in the appendix states that the family {PJk ^'^'y^Pt^^t'^^^TTt^^n > 2,k ^ 
{1, . . . ,n - l},m G {1, . . . , - l},s G [0, 1/n]) is bounded in ^|q,| +4(M'^). Since k < [nt\ 
when (t^, t) G 71, this imphes that {p'^A~ ,n > 2,k e {I, . . . , n — 1}) is 0{l/n) in +2(K'^)- 
Let us now prove the same for p"'"*". After \a\ integrations by parts and after setting 
z = y — C\/t — s, we get that — p) *g^a p)+ is a weighted sum of terms of the form 
^{Xa'ua2) - see Lemma [15]- with 

2 C) = ir+(5' - - CVt^s) 

dTiPn-p){s,x,y - CVt - s)p{t -s,y- (Vt - s,y) 

and |ai| + \a2\ = |a|. Now, from Corollary [22] in the appendix, {pn—p,n > 1) is 0(l/n) in 
G2{^'^) so that, using the same arguments as in Step 2 of the proof of Proposition [5]-(i), we 
get that {{pn -p) *g,aP)+ is 0(l/n) in H\a\ +20^'^) ■ Since p^A^{t,x,y) = 1^+ (t^, t)((p„ - 
p) *g,a p)+{tk,t,x,y), we conclude that {p^,l^,n > 2,k G {l,...,n — 1}) is 0{l/n) in 

Lastly, starting from (|43p , Lemmas [21] and [23] with j = 2 imply 

Lemma 20. Under (B) and (C), there exists a 0{\/v?) sequence (vTg,?! > 1) in ^4(M'^) 
such that for all t G (0, 1], / G C;;^i(M<^) and x G M^, 

Altf{x)= [ f{y)7Tnt,x,y)dy. 

Statement (iii) of Theorem [6] is now proved: it follows from (|45p . (|44p and Lemmas [TT] 
[m [IE and [20] 

We now also have all the tools to prove Theorem [8] Indeed, note that (|38p combined 
with Lemmas [21] and [23] imply that we have an expansion of arbitrary order j for pn — P- 

i=2 ' k=0 i=2 

Since {r"''^,n > 1) is 0{l/n^) in ^2?(]R'^) and {ip^u ,n > 1, G {0, . . . is bounded in 
^2i(K'^), this gives ([23|) with (7r„,i,n > 1) bounded in a2i-2(K'^) and > 1) and 

« ^,n > 1) bounded in g2i{^'^)- 
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4. Appendix 

4.1. Kernel of K^'^ . Here we make explicit the kernel of the remainder i?"'"', recall (|39p : 

Lemma 21. Under (B) and (C), for each j £ W, there exists a 0{l/n^) sequence 
(r"'^n > 1) in G2j(^'^) such that for all t G (0, 1], n > 1, / G C^^iR'^) and x G R'^, 



R?'f{x)= / f{y)r^^^{t,x,y)dy. 
Proof. From ([39]) and ([371), i?"'-'' = i?"'/ + R'lj where 

nt—\nt\/n f'Si rSj 

Let us first deal with R^'l- Using the fact that /c > 1 for the first equality, ([52]) for 
the second one, the fact that = Pt-t^-s for the third one, and ([Ml) and 

Fubini's theorem for the last one, we have for all / G C'^^i{R'^), x G M'^, t G [0, 1], n > 1, 
A: G {1, . . . , [nt\ - 1} and s G (0, 1/n), 

= [ Pn{tlx,Z,)'^y^^^Pt_t^ f{z,)dz, 

= ^ Yl IPnitl,x,zi)gj^o.Azi)P^{hj^^^id"Py^_,)Pt-ti^J{zi)dzi 

= yZ Pn{tl,X,Zi)gj^a,l{zi)Pn{s,Zi,Z2)hj^a,l{z2)0°Pt-q-sf{z2)dZ2dzi 
l<\a\<2j 1=1 -^^^ -^^^ 

f{y)Vq^{s,t,x,y) dy 
where c^"^-' = Ei<|a|<2j El='r V'^^a.i ^i^h 

ni,a,As^'t,X,y) = l(o_i)(s)l[in^^,i](i) 

(f'ri *gj,c,i,0l5n)(tfc,ife + S,X,Z2)hj^a,l{^2)d2P{t - - S,Z2,y) dz2. 



Now, setting q2l\,{u,x,z) = lu^ ,^{u){pn ^ ifi Pn){tl,u, x, z), it follows from Propo- 

sition [5]-(ii) that io^rl'^ai^'^ ^ E {l,...,n}) is a bounded sequence in ^(R^). Since 

'^t£a,ii^^''^^^^y) = l(o,^)(^)l[tI?+i,i](^)(Cf^a,/ + ■5'*'^'^)' Proposition [5]-(i) 

shows that {f^n^ijU > l,k G {l,...,n}) is bounded in 'H|q,|(M'^), so that {ip^n,n > 

1, G {1, . . . ,n}) is bounded in W2j(M'^)- 

When A; = 0, we have in the same way for all / G C|^i(M'^) 



'^7i/nPt-i/nfix) = / fiy)ipQ'\s,t,x,y)dy 



22 JULIEN GUYON 

where (^q'-' = Ei<\a\<2j Ya=i '^l^i with 

'/'O,!/!*' ^' y) = l(0,i)('5)l[i 1] it)9j,a,l{x){Pn *hj^^^,,a P){s, t, X, y). 

Again Proposition [5]-(i) imply that {ip^'\n > 1) is bounded in 7{2j{^'^)- 

Eventually, for all / G C^^ii^'^), we have Riif{x) = f{y)rl'\t,x,y) dy with 

ri'^{t,x,y)= y2 / / Vt^^^ i^j+i^''^^^^y) dsj+i-- -(132(181, 

f^^Q Jo Jo Jo 

and since the family "'^ ^, n > 1, A; G {0, . . . , n}) is bounded in T-L2j+2(^'^)-i the sequence 
(r"'^n > 1) is 0{l/ni) in C?2i(M'^)- 
As for -^2 * ' similar arguments lead to 

PVnt\/nKlvnt\/nfi^) = ^ f {y)^'' {s , t, X , y) dy 
H<2j 



where </."'^' = Ei<|a|<2,- ^1=^ with 



[nt\ [nt\ 



n n 



+ S,X,Z2 



hj,a,l{z2)d2P [t S,Z2,y \ dz2 + l{o<^<t< Ijfi'i.a,/ (a^) {pn *hj,^^i,aP) {s,t,X,y). 

We can treat (p'^'j exactly as we have treated 9?"n"'„ and get that {(p'^'^ ,n > 1) is bounded 
in 7i2j{M.'^), so that Rj^'^ has a kernel (r2'"',ra > 1) defined by 

r2'\t,x,y)= / ••• / (j)'''^+'^{sj+i,t,x,y) dsj+i---ds2dsi 

Jo Jo Jo 

which is 0{l/n^) in g2j{M.'^). 

Eventually, putting r"'-' = r"'-' + '"' completes the proof. □ 

In particular we have 

Corollary 22. Under (B) and (C), {pn-p,n> 1) is 0(l/n) in ^2(1^'^)- 

Proof. From (j38p applied with j = 1 and Lemma [211 we have for all / G C|^[(M'^) 

fiy)ipn - P)it,x,y) dy = A^f{x) = R'^'^fix) = [ f{y)r^^\t,x,y)dy 

jRd 

so that Pn — P = r"'^, and Lemma [2T] gives the result. □ 
Eventually, we have kernels for the operators P!lL*,Pt-t": 

h J k 

Lemma 23. Under (B) and (C), for each j G N*, there exists a hounded sequence 
('(/'"n'^n > 1, A; G {0, . . . ,n}) m^2j(K'^) such that for a// 1 G (0, 1], n > 1, A: G {0, . . . , [nt\}, 
f G C^iCM"^) and x G W^, 

P^L*Pt^tif{x)= [ fiym'{t,x,y)dy. 

jRd 

The proof is omitted since it copies the arguments of the proof of Lemma [21]- it is even 
a bit simpler. 
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4.2. Operators on Qi{W). When tt G ^K^K ), 7r(t, € L'' 
and n > 1 we can define two functions PsVr and P^vr on (0, 1] x 



so that for s G [0, 1] 
< IR^byP,^(t,-,y) = 



l{.<t}-Ps(vr(t,-,y)) and P^iT{t,-,y) = l{s<t}Psi-^it,-,y)), i-e. 

(51) P,7r(t,x,y) = l|,<t}E[^(t,Xj,y)] and PXt, y) = l|,<i|E [vr (t, y)] . 
We also write A^vr = P^vr - P^tt. For j G N* we denote by the family , ,n > 



l,s E [0, 1/n]) of operators on GiO^ ) defined as in ([3T]) by 

<i>^;,^/„7r(t, X, y) = E [L|Pi/„_,vr {t, X^'\ y) 

i.e., using ([30]) . 



(52) 



s,l/n 



l<|a|<2i 1=1 



Denoting by Cb{Gi{W^),gv{W^)) the space of all morphisms mapping any bounded subset 
of Qi{W^) into a bounded subset of QuiW^), we then have 

Proposition 24. Under (B) and (C), (P^,s G [0,1]) and (PJ',s G [0,1], n > 1) are 
bounded families in Cb{Qi{'K'^)), and is a bounded family in ^^(^((M'^), g,+2i(IK'^))- 



Proof. Let us first deal with (P^). Let vr G ^;(M ). P^ is measurable. Moreover, Lebesgue's 
dominated convergence theorem shows that Ps7r{t,x, •) is infinitely differentiable and that 
for all /3 G N'^ 



d^PMt,x,y) = l 



{s<t} 



E 



Hypothesis (A) ensures that a version of can be chosen such that for each t > 0, the map 
X I— > X^ is infinitely differentiable (see, for example, [lO]). Since d^7r{t,-,y) G C^i(M'^), 
it follows from Theorem 3.14 page 16 in [H] that dyPs7r{t, -,!/) is infinitely differentiable 
and that for all a G N'^ there exists universal polynomials (Hq,^^, < |a|) such that 

(53) d^d^PMt, X, y) = l{,<t} E ^ [^2 5f ^ (i> X^s^v) na,^ (^^^^ l^^l < !«!) 

H<l"l 

with 

(54) sup E[n„,^ (5^X,^ \v\ < \a\f] < oo 

se[0,l],i'GK'' 

for all [//[ < |a|. As a consequence, Ps7r(t, •, •) is infinitely differentiable and using Cauchy- 
Schwarz's inequality, (jl4p and (j54p . we see that for all bounded B C aK^'') and a,(3 G N'^, 
there exists two constants ci > and C2 > such that for all vr G ;S, s G [0, 1], t G (0, 1] 

and X, y G M'^, 



(55) 



d^d^Ps7:{t,x,y) < cil|,<i}t-(l"l + l^l+'^+')/2E 



exp 



-C2||X^-2/||Vt 



1/2 



Now, partitioning Q into — y|| < \\x — y\\ /2} and {\\X^ — y\\ > \\x — y\\ /2}, we have 



(56) E [exp(-C2||X|-y||Vi) 



< 



y\\ < - y|| /2) +exp ( -C2 ||x - y|| /At 
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Using (fT6|) for p G Q{W^) for the fourth inequahty, we can find 03,05 > and 04,05 > 
such that for all s G (0, 1] and x, y G W^, 

¥{\\X^, -y\\<\\x-y\\/2) < F {\\X^ - x\\ > \\x - y\\ /2) 

= / ^{\\z-x\\>\\x-y\\/2}Pis,X,z) dz 

= / l{||^||>||^_j^||/2^}P(s,X,X + C\/s)s'^/^ 

- ^3 / l{||5||>||x-j;||/2Vi}exp(-04||Cf )dC 

(57) < osexp ( -06 ||x - /s 



Eventually, from (156]) and (f57l) . we can find 07 > and cg > such that for all s £ [0, 1], 
t G (0, 1] and x,y £ R^, 



l{s<t}E 



exp (^-02 - /t^ < C5exp (^-06 ||x - /t^ + exp (^-02 ||x - /4tj 
(58) < C7exp ( -08 ||x - y||^ /t 



It is enough to inject ([58|) into ([55]) to complete the proof for (Ps)- 

This proof naturally extends to the case of (-P"). Indeed, ([53]) holds with (X",P") 
instead of {X,P). Moreover, from Lemma [26l ()54p holds uniformly in n with instead 
of X. Eventually, (j57p holds with X"' instead of X, uniformly in n because {pn,n > 1) is 
bounded in g{R'^). 

As for , it is enough to use (j52p . the boundedness of (P^) and (-P"), Remark 1111 and 
the facts that multiplication by a function in B belongs to Ci,{Qi(M.'^),Qi{M.'^)) and that 



4.3. Moments for the Euler scheme and its derivatives. Let us assume (A). Then 
it is known that X"'^ has bounded moments of any order and that for all g G N, one can 
find > such that for all x G M"^, 

(59) sup E [IIX"'^!'^] < c(l + ||xf ) 

te[o,i],n>i 

(see 115]). From x ^ X^'^ is infinitely differentiable and we shall see that analogous 
upper bounds hold for its derivatives. Following [11], for m > 1, we denote by xj:™''^''^'^ 
the m-th derivative of x 1-^ X^'^ at point x. It should be thought of as a d x matrix. 
For instance, X^^^'"''^ is the jacobian matrix of x 1-^ X"'^. Differentiating Q, we have 



(60) 



/■* (1)^ )xl')'^fds + y t af\x'^'\, )x\''^f;dBl 
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where I stands for the identity matrix and aj is the j-th column of a. Besides, by induction, 
there are for each m > 2 universal polynomials Pm,,j, j £ {0, . . . , r}, such that 



(61) X, 



{m),n,x 



+ r<3"'ni 1/ ds+y Tq"'", ds^', 

^m,0,lns\/n / j ^■tn,],\ns\/n si 



where 



(62) 



'm,j,t - Prn,j[Crj (,Aj (,A^ j,Aj , • • • , Aj j. 

This is analogous to (1.8) page 4 in pjj. Then we have 

Lemma 25. Under (A), for all m > 1 and q G'N, there exists c > and g' G N such that 
for all X £ R"^, 



(63) 



sup E 

iG[0,l],n>l 



{m),n,x 



< C 1 + llxl 



Proof. We give a proof by induction on m. Let us first assume that m = 1. Let q G N. 
Prom (j60p . and observing that (A) states that b^^^ and all the trj^^ are bounded, Jensen's 
and Burkholder-Davis-Gundy's inequalities lead to the existence of c > such that for all 

t e [0, 1], n > 1 and x € M"', 

rt 



E 



{l),n,x 



< C 1 + 



E 







X 



{l),n,x 
\ns\/n 



ds 



Taking this inequality at time \nt\ jn and applying Gronwall's lemma, we get that 



sup 

te[0,l],n>l,a: 



E 



X 



(l),n,a: 
[nt J / n 



< OO. 



From ([U, one easily checks that the same holds at time t instead of [nt\/n, so that 
holds for m = 1 with q' = 0. 

Let us now assume that ()63p holds for the m — 1 first derivatives. Let q £ N. From 
([6T|) . and observing again that (A) states that b'^^^ and all the crj^^ are bounded, Jensen's 
and Burkholder-Davis-Gundy's inequalities lead to the existence of ci > such that for 
ah t G [0, 1], n > 1 and x G W^, 



(64) E 



X 



{m),n,x 



<Cl 




X 



(m),n,x 
[nsj /n 



•^0 j=0 



ds 



Using (j62|) . the induction hypothesis, (A) and (|59|) . we find C2 > and q' £ N such that 
for aU s G [0, 1], n > 1 and x G M*^, 



r 

j=0 



m j', [nsj /n 



< C2 1 + ||a;| 



Thus, taking (j64p at time [ntj/n and applying Gronwall's lemma, we find c > such that 
for ah X G M*^, 



sup E 

tg[0,l],?i>l 



X 



(m),n,x 
[ntj /n 



< C 1 llxl 
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From one easily checks that the same holds at time t instead of [nt\ /n, which completes 
the proof. □ 

Observe that, under (B), the above proof holds with q' = so that we have 

Lemma 26. Under (B), for all m> 1 and g S N, 



sup E 

te[0,l],n>l,x( 



< OO. 
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